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5^ . Abstract 

In a dilute system of instantons and antiinstantons, the Ua(1) and scale anomalies are 
• 1—1 , 

^ ■ shown to be directly related to the bulk susceptibility and compressibility of the system. 

H ■ 

Using 1/Nc (where Nc is the number of colors) as a book-keeping argument, mesonic, 
baryonic and gluonic correlators are worked out in p-space and Fourier transformed to 
x-space for a comparison with recently simulated correlators. The results are in overall 
agreement with simulations and lattice calculations, for distances up to 1.5 fm, despite 
the fact that some channels lack the necessary physical singularities. We analyze various 
space-like form factors of the nucleon and show that they are amenable to constituent 
quark form factors to leading order in l/Nc- Issues related to the lack of confinement in 
the model and its consequence on the various correlation functions and form factors are 
also discussed. 
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1 Introduction. 



An outstanding problem in QCD is the understanding of the hadronic spectrum from 
first principles. Decades of dedicated lattice simulations have shown that the problem 
is difficult when all QCD degrees of freedom are taken into account. Through the years 
there have been numerous proposals, both theoretical and numerical, which suggest that 
only some relevant degrees of freedom may be important for the bulk aspects of the 
hadronic spectrum. The proposals range from lattice cooling procedures to semi- 
classical techniques [0]. 

Recent lattice simulations based on cooling procedures have suggested that instanton 
and antiinstanton configurations may account for a large part of the hadronic correlations 
0, ^, although the local character of the cooling algorithms may not totally rule out 
persistent quantum effects at large distances p. On a periodic lattice (without twists), 
instanton or antiinstanton configurations are necessarily singular [Q. Their continuum 
analogs are the BPST instantons in singular gauge, a point of some recent concern 0. 
Notwithstanding such concerns, an impressive amount of results, both from cooled lattice 
simulations f^, Q and from random instanton simulations seem to indicate that the 
basic features of the hadronic spectrum may emerge from a dilute ensemble of singular 
instantons and antiinstantons. 

Sometime ago, 't Hooft's suggested that instantons provide the answer to the axial 
Ua(1) problem. In the presence of instantons or antiinstantons, light quarks acquire 
zero modes, which bunch into fiavour-singlet configurations ('t Hooft's vertices) thereby 
dynamically breaking the axial Ua(1) symmetry. At low energy, 't Hooft's interactions 
provide interesting correlations in various hadronic channels, as noted by Callan, Dashen 
and Gross 0|, and analyzed using QCD sum rules |Ty], resummation procedures [|1^, [13 



instanton simulations [T^], and bosonisation techniques |15, 16, [T^) HH- 

In this paper, we will assume that the QCD partition function simplifies into a grand 
canonical ensemble composed of 't Hooft's vertices, with an apriorily unspecified measure 
for the instanton-antiinstanton interactions. We will further assume that the ensemble is 
dilute with a screened topological charge, as discussed in Refs. |TB|, [T^, ^ . The screening 
is expected from the feedback of the light quarks on the instantons and antiinstantons 



in the vacuum In this respect, problems related to the original choice of the 



instanton-antiinstanton ansatz 12311, as well as the limitations associated with the stream- 



line approach ||2^, are somehow irrelevant. Analytical and numerical calculations with 



such an ensemble have led to a satisfactory phenomenology p, |T3|, [19 1 



The purpose of this paper is to show that the analytically derived results using either 
resummation techniques [0 or bosonisation techniques ||15|, |19| for two and three flavours 



in momentum space, are consistent with the recent four-dimensional simulations as well 
as with the cooled and quenched lattice simulations up to a distance of 1.5 fm. At 
larger distances, the lack of confinement shows up in the form of spurious oscillations. 
The physics of a screened gas of instantons and antiinstantons is well described by simple 
mean field arguments |]I5|, In section 2, we recall the effective action for a random 



instanton gas using an approximate bosonisation scheme. In section 3, we discuss the 
structure of the massive quark propagator both in momentum and coordinate spaces, and 
comment on heavy-light correlators. We note that in the long wavelength limit, the quark 
propagator becomes tachyonic for all current quark masses. In section 4, we outline the 
result for the quark condensate in the random instanton gas. In section 5, we give a 
brief account of the various mesonic correlators, including the scalars. We discuss issues 
related to the mixing between the scalars and the fluctuations in the instanton scalar 
density through the scale anomaly. The mixing between the pseudoscalar singlet and the 
fluctuations in the instanton pseudoscalar density yields naturally to a resolution of the 
Ua(1) problem. Issues related to the 77 — 77' mixing are also discussed. The p-space results 
are discussed in detail for scalars, pseudoscalars and vectors, and compared with the x- 
space simulations. While the analysis of the p-space pseudoscalar correlator shows clear 
evidence of poles, the vector correlators simply exhibit two constituent quarks. In section 
6, we briefly discuss nonstrange baryons in x-space. In section 7, scalar and pseudoscalar 
gluonic correlation functions are discussed. In section 8, quark and gluon form factors of 
the constituent quark are discussed. To leading order in l/A^^c, they saturate the nucleon 
form factor following from the point-to-point correlator in terms of loffe's current. Our 
conclusions and recommendations are summarised in section 9. 

The details of the bosonisation techniques are given in Appendix A. In Appendix B, 
we provide a direct calculation of the quark condensate. In Appendix C, the necessary 
elements for a Gaussian approximation are presented. In Appendix D, the effective action 



for the singlet and octet pseudoscalars is explicitly worked out. In Appendix E, an exten- 
sion bosonization scheme is presented. In Appendix F, the various expressions entering 
the unconnected parts of the mesonic correlators are summarized. Some of the difficul- 
ties related with the expansion of the mesonic vertices involving the strange quark mass 
are discussed in Appendix G. In Appendix H, we outline the essentials of our numerical 
procedures. 



2 Model 

• Effective action 

t'Hooft has shown that at scales larger than a typical instanton size p (fixed throughout 
this paper), instantons induce fiavour mixing between the light u,d and s quarks in the 
form of determinantal interactions (t'Hooft determinants) ^ 

det± = ^det;, (m^V - |(/ V^^^^V^ / <P^^S^'r)) (1) 

where m-^^ = diag(m, m,ms) is the current mass matrix for {u,d,s) quarks, (f)± are the 
instanton-antinstanton zero modes, ip is the fermion field in the long wavelength limit, 
and 5*0^^ = — (i^ + im) is the free fermion propagator. The averaging implied by (■ ■ ■) is 
over the instanton and antiinstanton color orientations. 

A random system of instantons and antiinstantons that is compatible with the Ua(1) 
and scale anomaly yields the generating functional ||, ^7\, 

Z[7],r]^] = J dn+dn_V^V^^ /u(n+,n_) e-/^'^^['''''^'"+'"-l (2) 

where 

C[ri, ri\ ?T,_|_, n_] = if)'^ Sq'^iIj — log det+ — n_ log det_ — ^/'^r/ — r/V (3) 

in (^ at the saddle points. Throughout this paper, the generating functional will be used 
to carry calculations to leading order in l/N^, where A^^^ counts the number of colors. The 
counting will be understood just as a convenient way of organizing the calculation, with 
Nc = 3. In the presence of instantons, conventional Nc arguments have to be amended 
25|1 (see also below). 



We are using a coarse grained action for the description of the instantons and anti- 
instantons, as discussed in [^. The relation to the uncoarse-grained approach, follows 
from the identification 

N± 

n^iz) = ±J26\z-z,) (4) 

i=l 

at the saddle points {N^ = 3 >> 1). The coarse grained version highlights the role of the 
scalar and pseudoscalar glueball fields, and their mixing to the quark-antiquark excita- 
tions. The measure /i(n+,n_) refers to the distribution of instantons and antiinstantons 
in the vacuum without the quarks (quenched approximation). Its form is generic, and 



follows solely from the Ua(1) and scale anomalies pO 



/i(n+,n_) = exp[^ - — J d z{n^{z) + n {z))\^og - 

- ^Jd'z{n^z)-n-{z)f^ (5) 
where n = N/V4 ~ N^. is the mean instanton and antiinstanton density in the thermody- 



namical limit, x* = ^ l23[ the quenched topological susceptibility 

2 



X*=< (^J d^z{n+ -n^){z)^ >Nf=o (6) 



and (T* the quenched compressibility (a = 1, 2, ...) 

r2\ a-1 



=^<[jd'z{n^ + n--n) {z)) \^^=o (7) 



Anomalies 



The coarse grained effective action along with the measure @ satisfies both the 
axial U(l) and scale anomaly ^0\. Indeed, in the chiral limit the determinants in 



acquire a phase under a U(l) axial rotation, hence a non-conserved axial-singlet current, 

d^j^,iz) = 2Nf{n+ - n_){z) + 2i Tr^ m 'iP\z)^^'iP{z) (8) 

Also, in the quenched approximation, the measure (j^) is not scale invariant. As a result, 
the divergence of the dilatational current (trace of the energy momentum tensor 0^,^), is 
not conserved. 

An 2 
e,,(z) = (n+(z) + n-{z)) + -(n+(z) - n.{z) f + 0{Ns) (9) 



Comparison with the QCD form of the trace anomaly |26| gives crl/n 12/llNc ~ l/A^^c. 



In the fermionic part, we note that for n+ = n_, the generating functional @ involves 
only the combination (det+det_) and is invariant under ?7l(3) x f/ij(3). The fluctuations 
in (n+ + n~) involve mixing between the isosinglet scalar and the scalar "glueballs", 
through 

- ^ y c/S(n+ + n-){z) ln(^det+ x det_^ (10) 

For ^ n_, the combination (det+/det_) is also allowed. The latter dynamically breaks 
the axial U(l) symmetry through 

- i + n-)Wln(|^±) (11) 

as originally suggested by t'Hooft. The fluctuations in {n^ — n^) will mix with the 
isosinglet pseudoscalar, thereby resolving the Ua(1) problem (see below). 

• Bosonization 

In vacuum, the packing fraction is given by the dimensionless combination 

„y._^,*.10- (12) 

Since the density n ~ Nc, the packing fraction is of order A^^. The value (|1^) is small, 
and an expansion in the density is justified except in the presence of infrared singularities, 
as we will specify below. In this spirit, the generating functional (0) can be bosonized 
approximately by inserting the identity [] 

1 = Jv7r^VP^exp(TTf Jdkdl P^{k,l)(n^{k,l) -e^{k,l)^^ (13) 

in the partition function (|^), where vr^ and P^ are bilocal auxiliary fields and Nf x Nf 
valued such that P^{k, I) = P^{k — I) and similarly for 7r^(fc, /). Also (see appendix A) 

9^{k, I) = {^\k)S^'<p^{k)<j>^\l)S,'m) (14) 



The trace over flavor indices is understood in the exponent (|T3|) . The auxiliary fields vr^ 
can be eliminated by using the saddle point approximation. From Appendix A, we have 

Z[r^, ryt] = J pp± e~r,ts[p+,P-]v e-Scff(P±) (15) 



^In what follows, we use the shorthand notation dk = d^k/{2Ti)'^ and dz = d'^z when integrating. 



where the effective action is given by 

4 „, _ „ , 4 



Ti f / A 4 

N,Ti(\ogS-^\P+,p-]) + - dz [Ti \n—P+(z)+TT \n—P-(z) 

2 J \ np np 

2 J dzTifm{P+{z) + P'{z)) (16) 



The trace Tr is over ffavour and Dirac indices as well as four momenta, the Tr^ over flavour 
indices and the det is over flavour indices as well position space. We have explicitly used 
n+ = n_ = n/2 and defined the operator momentum dependent inverse propagator 

S-'[P\ p-] = f-tm- tJm + P-75i fl - ^^fl \A4 (17) 



Here, k refers to the momentum operator and is the induced momentum dependent 
screening quark mass. The screening mass arises from averaging over the instanton zero 
modes [1^. From Appendix C, we have after rescaling 



Mk = Mkim) = A(m) — A;V2 = A(0)— A;^ ^ V_ (l,Ko - IiK,)j (18) 



where t = kp/2 is the argument of the McDonald functions / and K. From (|T6]) and ([T7|) 
it follows that in the long wavelength limit quarks in the instanton vacuum interact via 
the exchange of effective bosonic fields P^. The latter are Nf x Nf valued and may be 
parametrized as 

Other parametrizations are also possible [^]. We note, however, that to the order we 
will discuss the correlation functions below (gaussian approximation), the results are 
parametrization independent. 

• Gap Equation 

The matrices k and a are Nf x Nf valued and hermitean. The k variables can be 
identified as the pseudoscalar Goldstone modes, except for kq. The matrix a contains 
the massive scalar-isoscalar and scalar-isovector excitations. The non zero value of a in 



^This constraint will be relaxed below to address the rj' mass and the gluon correlators. 
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the vacuum follows from the saddle point approximation to (|T6D by setting k = 0, and 
switching off the sources. The result is an integral equation for each flavour 



^ ' n J (27r)4 - 2mMk{m)k^ + (P + m^)Mi{m) ^ ' 



where we have set Mfc(m)/Mfc(0) = A(m)/A(0). For an instanton density n = 1 fm~^ 
and size p = 0.33 fm, the behaviour of the constituent quark mass Mk{m) versus the 
dimensionless combination z = kp/ 2 is shown in Fig. 1 with current quark masses m = 
(dashed), 5 MeV (dotted) and 10 MeV (solid), respectively. For k ^ 0, Mo(m) — >■ 
A (m) (n/A'^c)(2vrp)^, while for k ^ 1, Mk{m) falls off like 1/k^. The width at half maximum 
is of order 1/p. The result ( pO]) was also obtained in using different arguments. 



3 Quark Propagator 

In a random instanton gas, quarks are "screened" . The light fermion propagator acquires 
a momentum dependent mass 0. For one flavour, the results of Appendix A give (unless 
specified, we denote Mfc(m) by M^) 



S{k,m)-^_^^ ^_,(^_p/M,) ^^^^ 



In the massless case 13 



We note, however, that at low momentum, Mk{m) acquires a non-analytical contribution 
(k^O) 



Mkim) ~ Mo(m) (l + 3^^ log -e'+^'j (23) 

V 2 / z=kp/2 



^In general, the quark propagator is gauge dependent. Our case is no exception, and the present 
discussion should be understood as the evaluation of the quark propagator in a random and classical 
background of instantons and antiinstantons in a singular gauge. 
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where c = 0.577 is Euler's constant. The screening mass Mo(m) does not show up as a 
simple pole. What this means is that, as A; — > 0, the screened quarks become tachyonic. To 
the extent that long wavelength quarks are unphysical, this should be of no real concern. 
However, since the instanton model does not provide for confinement, these "unphysical" 
effects will contaminate all large distance behaviours. The instanton simulations 0, or 
cooled lattice calculations have not probed large distances. 

With this in mind, we now proceed to x-space with the decomposition 



S{x, m) = Sq{x, m) + S'i(x, m) 



(24) 



where S'(x, m) is the Fourier transform of ( pT)) and 



So{x,m) 



im 
4^ 



x^ 



K2(mx) H — Kiimx) 



X 



(25) 



is the free propagator of quark of mass m. 5*1 (x,m) follows from (|2^) and will be under- 
stood as 



Si{x,m) =i(i:Sf'^-ST 



(26) 



where 



5p(a;, m) 



1 d 



+ 00 



dk 



X dx \A'k'^x Jo k"^ + {m — k'^/M^) 



;Ji{kx) 



(27) 



and 



ST^'^ix.m) 



1 /■+-,, ^(m-fcVMfc) 

dk— -, . „ , Ji(fcx) 



A'k'^x Jo k"^ + {m- k'^/Mk)' 



(2^ 



Figs. 2a and 2b show the behaviour of mTrS'(a;, m)/TrS'o(x, m = 0) (chirality flip) and 
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Tr74>S'(x, m)/Tr745'o(a;, m = 0) (chirality non-flip) versus x up to 2 fm, for quark masses 
of 5 MeV (lower curve) and 10 MeV (upper curve), respectively. The squares refer to the 
results of simulations of Ref. |§] using 128 instantons and 128 antiinstantons in a periodic 
Euclidean box of 3.36^ x 6.72 fm^. These simulations were carried out with equal u and 
d quark masses of m = 10 MeV. The small discrepancy in the chirality flip part of the 
propagator may be due to the fact that the instanton simulations make explicit use of 
the single-instanton distorted propagator for the nonzero mode part, while the bosonized 
constructions presented above make use of the undistorted light quark propagator for the 
nonzero mode part. 

Figs. 2c and 2d show the chirality flip and non-flip part of the quark propagator 
for m„ = 10 MeV and rris = 140 MeV over a wider range of x. The larger the quark 
mass, the larger the oscillation in the quark propagator at large distances. These spurious 
oscillations are due to the appearance of the tachyonic mass (|23| ) and the occurence of 
the combination (m — k'^/Mk) in the quark propagator, and will cause most correlators 
to lack scaling at large distances (typically of the order of 2.5 fm and larger) as we will 
discuss below. We have checked that these oscillations persist in the massless case. In 
fact for m = Fig. 2c is almost unchanged. 

At this stage, we should point out that our treatment of the current masses is only 
approximate, given our definitions ([^) and (§). We will check below that the linear effects 
in the current mass do reproduce known results, while the non-linear effects cancel out 
at large distances, leaving us with the expected masses for the strange pseudoscalars. 
Similar observations apply to the instanton simulations in |^, although the handling of 
the current masses is not necessarily the same as the one discussed here. 

Finally, we note that the naive interpretation that the x-space version of Tr(5'(x, m) (1± 
74) /2) as the correlator of a light quark in the field of an infinitely heavy quark p9[] 
overlooks the issue of binding. As it stands, the non-relativistic projection of the heavy- 
light-propagator without the Wilson line (Coulomb field) for the heavy particle reflects 
solely on a screened light quark. In a heavy-light system like a D or B meson, the light 
quark is expected to bind to the heavy source, causing the spectral function to develop a 
pole instead of a cut. A detailed analysis of systems with few heavy and light quarks in a 
random instanton gas has been given in Ref. |^ . In the Coulomb field of a heavy quark, 

11 



the light quarks bind with a binding energy of the order of a quarter of the screening mass 



4 Quark Condensate 

The formation of a quark condensate in the instanton vacuum follows from the random 
nature of the system. From our bosonized construction, the quark condensate is obtained 
from the effective action (0) through 

Since the present treatment is semi-classical, all the ambiguities associated with the cur- 
rent mass singularities are ignored. At the saddle point, a straightforward calculation in 
the m ^ limit gives 



Using the mass gap equation for zero current mass, the term in brackets multiplying A'(0) 
vanishes and we are left with 



< >= -AN, { - I ^^^^^^^-J^ (31) 



^2N, ^ '~ J (27r)4 (P + Ml 

As a check, we show in Appendix B how this result can be recovered from the origi- 
nal definition in the saddle point approximation, prior to the bosonisation procedure. 
Numerically^, 

^(0) = 2 / IK-T.TTJTI^TK (32) 



IN, ' J (27r)4(F + M|) 

so that 

d^k M, 

which is the expected result to leading order in l/N,. 



< >= -AN, I , = -(5(0, m 0+)) (33) 



*With our choice of parameters, the discrepancy is (10) MeV^. 
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5 Mesonic Correlators 



To leading order in l/Nc-, the mesonic correlation functions follow from (|T6|) by differen- 
tiation with respect to the external sources in the presence of the auxiliary bosonic fields 
P^. Generically, 

C^{x) = {T*ij+-fij{x) V^T^/'IO)) (34) 

with 7 = (1,75,7^,757^''^/^,'^) ® i^^T"). From (0) we have = C° + C^, where the 
connected part of the correlator is given by 

C',{x) = -^^|l)P±Tr(S[x,0;P]7S[0,x;P]7)e-^-[^^l (35) 

and the unconnected part is given by 

C'.ix) = ^ji^|pP±Tr(S[a;,a;;P]7)Tr(S[0,0;P]7)e-^-[^^l (36) 

Typical diagrams contributing to (|3^) are shown in Figs. 3a and 3b. Only the diagrams 
in Figs. 3c and 3d are dominant. They will be the only ones discussed here. In p-space, 
the contribution of Fig. 3c is 



C°(P) = -Nc I ^Tr (5(1)75(2)7) (37) 



where 5(1,2) = S{k ±p/2,mi.2) for two arbitrary flavors. This contribution in the long 
wavelength limit reflects on the lack of confinement in the model. The contribution of 
Fig. 3d is 
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The extraction of A± and R± from ([l6|) is performed in Appendices B and C, respectively. 
With the above approximation in mind, the total correlation function reduces to the sum 
of (|37[) and (|38D , i.e. = C° + C;J. The results ( ^7| - ^ ) were first derived by Dyakonov 
and Petrov for two massless flavours using detailed resummation procedures . 



At this stage, it is interesting to compare the expression we have for the mesonic 
correlator in the instanton model with the one derived in planar QCD2. In the large Nc 
limit, the two-fermion cut in QCD2 is infrared sensitive and cancels exactly against the 
infrared sensitive one-gluon exchange graph This cancellation makes explicit use of 
Ward identities in Feynman graphs. It is essentially quantum and thus absent from the 
present semiclassical argument. The lack of confinement in our case will have dramatic 
consequences on the large distance behavior of the various correlation functions as we will 
discuss below. 

The expressions used to generate the various correlators in p-space are tabulated in 
Appendix D. In Figs. 4, we show the behaviour of the connected (minus the vacuum) 
correlators in the various channels versus the momentum p, for the Axial- (A), Vector- 
(V), Pseudoscalar- (P), Scalar- (S) and Tensor-channel (T), without strangeness (Fig. 
4a) and with strangeness (Fig. 4b). Similar correlators are shown in Figs. 5 for the 
unconnected part. By about p ~ 10 fm~^ the correlations are totally washed out. The 
plots are for u and d quark masses of 10 MeV and a strange quark mass of 140 MeV. 
Although unconventional, this choice of the current masses allow for a comparison with 
the numerical simulations of Ref. M. 



5.1 Gell-Mann-Oakes-Renner Relation 

In the pion channel, a pole is produced by the unconnected part of the correlator that 
lies well below the two-constituent quark cut. This is a good example of an infrared 
sensitive channel, where a simple expansion in the instanton density fails. The presence 
of small denominators through zero modes forces the resummation of an infinite string of 
terms of increasing powers in the instanton density, causing the correlation function to 
develop a pole. Using the small momentum expansion (see Appendix F), we have in the 
pseudoscalar channel 

14 



= J^i^~+P' + ^(^'^P')) (39) 



4N, 



to leading order in the current quark mass m. Above, the decay constant / satisfies 



(27r)4 (P + MlY 



and the pseudoscalar mass M_ is given by 



-M! = 2m -— A(m) - / -—— ^ (41) 



AN^ - \2N^ ' ' J (27r)4/c2 + M2 

To this order, the quark condensate is current mass independent and is given by (|3TD . 
Thus, 

f^ml = -2m <^ij> (42) 

which is the current algebra result derived by Gell-Mann, Oakes and Renner (GOR) 
P^ . For equal u and d quark masses with m = 5 MeV, we obtain m.^ = 158 MeV, 
< ipip >= —(249 MeV)^ and to leading order in the current mass /,r = 88 MeV. Similar 
results can be derived for K and 1], although the small momentum expansion is no longer 
valid for the unconnected part of the correlation function with a large strange quark mass. 
This point is further discussed in Appendix F. 

In the expansion discussed above, the consistency of the GOR result can be further 
checked by noting that the unconnected part in the pion channel reads 



(43) 



where 



The term in brackets in the expression for can be identified with the usual pseudoscalar 
strength (7^. From (|3ll - |33D , it follows that Qt, ^ 2 < ipip > j fT,. 
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5.2 Pseudoscalars 



• TT and K 

Figs. 6 and 7 show the behaviour of the pion and kaon correlators versus x, respec- 
tively, as they follow from (|37D and ( ^) by Fourier transforms. The upper curve is for 
m = 5 MeV, while the lower curve is for m = 10 MeV. The squares are the results of sim- 
ulations using 128 instantons and 128 antiinstantons in a (3.36^ x 6.72) fm ^ periodic box. 
The dotted circles are the results from cooled and quenched lattice gauge calculations on 
a 16^ X 24 lattice with a physical lattice spacing of 0.17 fm. 

The momentum dependent parts display a low-lying spurious cut at about 627 MeV, 
as well as a pole in the scalar and pseudoscalar channels. In Fig. 8, we display these two 
separate contributions to the pion channel for an average quark mass of 10 MeV. 

The asymptotic form of the correlation function in x-space is strongly influenced by 
the position of the pole in most channels. Indeed, the large distance behaviour produced 
by the pole is approximately of the form 



-M±x 



Cl{x ^ oo) ~ ^^^—^ ^ (45) 

4 (27rM±x)'/' 

which is to be compared with the contribution of two "regular" (not tachyonic) screened 
quarks: 



with Mq = Mo(m). Beyond 2 fm, the running mass {m — k'^/Mk) causes (|46|) to oscillate as 
shown in Figs. 9 and 10. These oscillations, however, are overpowered in the pion channel 
given the very large signal caused by the pion pole compared to the spurious cut (about 
100 : 1). We note that for m = 10 MeV, the screening mass for the two screened quarks is 
about 627 MeV. In (^), the pseudoscalar mass squared follows from the GOR relation 



M! = -(mi + ma) ^ ^ (47) 
16 



In the case of the pion, we plot in Fig. 11 the total correlator times x^l'^ for quark masses of 
5 and 10 MeV, respectively. The pion mass sets in at about 2.5 fm. From the asymptote, 
we read a slope of m.,^ ~ 157 and 215 MeV, respectively. The agreement of the slopes with 
the GOR result provides consistency checks on the various Fourier transforms performed. 
We stress that to read the masses through slopes requires a proper identification of the 
preexponent power (here A raw plot of the total correlator versus x does not show 

any scaling up to 10 fm ! 

A similar analysis for the kaon channel is shown in Fig. 12, where only the rescaled 
and unconnected part C^{x) is shown. The connected part oscillates at distances of the 
order of 2.5 fm and larger, as shown in Fig. 9b, for two strange quarks. In contrast to 
the pion channel, the ratio of the connected to unconnected parts in this case is about 
5:1. The linear fall off in Fig. 12 sets in between 2 and 3 fm. From the slope, we read 
mK = 490 MeV, for m = 5 MeV and vris — 140 MeV. We note that all the non-linearities 
in the strange quark mass cancel out to give a kaon mass that is compatible with the 
mass obtained by a naive use of the GOR relation, as indicated above. 

• 7] and 7]' 

In the r] and r]' channels, the situation is a bit more subtle because of mixing and the 
anomaly. First, let us follow the nonet decomposition used in Appendix D, for the singlet 
(ko) and the octet (^g) excitations. The connected part of the correlator in the (00,08,88) 
channels reads 

C\p) = - j dkTr (75A^-S(A;i, m)75A^5(A;2, m)) (48) 
with any of the singlet or octet U (3) generator. Specifically, 

C°(x)= ^C°(x,m„) + ^C°(x,m,) 

C'osix) = ^ (C°(x, mj - C^x, m,)) (49) 
with C°(x, rriu) and C°(a;, m^) the correlators of two screened uu and ss quarks, with 
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rriu = rrid = 10 MeV and = 140 MeV. The behaviour of (^9[) is shown in Fig. 13 
versus x. The oscillations seen in all channels beyond 2.5 fm are due to the spurious 
quark modes. 

The unconnected part of the correlators in the singlet, octet and mixed channels follow 
from the results of Appendix D. Since the tjq and 1]^ fields are integrated over, we can 
desentangle them by a unitary rotation of angle (A = A_ {k, m, m), = A_ {k, nig, rus) 
and z = x*Nf/Nc) 

where 



z . f, . . z^ 2z , . ..V^ 



X±{k) = A + A, + - ± (^(A - A, + + y (A, - A) j (51) 
at the expense of rotating the vertices (sources) as well. The result is 



f {Rocos9 + Rossm9f ^ (-i?o sing + i^pg cos ^)^ ^ 

- I x^ik) + X4k) J 



1 /(i?08cosg + i?8sing)2 (-i?o8sing + i?8cosg)2\ 



1 _ , . (-Ro cos 9 + Ro8 sin 9) (i?08 cos 6^ + i?8 sin 9) 
CM - 2iV, 

^ (--Rq sin 9 + -Rqs cos 9){-Rqs sin 6^ + -Rs cos 6*) 
+ ' X_ik) ^ ' 



The poles in the unconnected parts are just the rj and rj' masses, since we have rewritten 
the singlet and octet correlators in the rj and t]' basis |. At low momentum 9 ~ —13.1°, 

^ Since the diagonalization is momentum dependent, it is not possible to devise a local source that 
would trigger precisely the rj or rj' quantum numbers without mixing. 
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which is to be compared with 9 ~ —11.5° in [|T^. From the effective action of Appendix 
D, we conclude that to leading order in the current masses 

m + TUs 



X±{k ~ 0) 



± 



< ihih > -\ — 
2N, 2 

.rUs — m — , z.r, 2zms—m 

^ 2Nr 2' 3 2Nr 



1/2 



Then, 



fmj, =2iV, A+(A;~0) 



and 



r2 2 



2N^ X_{k ~ 0) 



(55) 



(56) 



(57) 



The above relations give m^/ = 1163 MeV and m.^ = 557 MeV. These values are to be 
compared with m^/ = 1172 MeV and = 527 MeV for < ipip >= (—255 MeV)^ and 
/ = 91 MeV as used in Ref. 0. From (M), ([ 



and (|57D we have 



f{ml,+ml-2mj,)=2NfX. 



(58) 



which is the Veneziano-Witten formula 25, pSl. 



Figs. 14a and 14b show the behaviour of the rescaled and unconnected parts of the 
correlators versus x. At about 3 fm, the asymptotic slopes set in. The large distance 
behaviour being dictated by the smaller pole, we obtain a slope of 220 MeV (essentially 
the pion mass) when the term x*Nf/Nc is switched off, and a slope of 480 MeV (essentially 
the 7] mass) when it is not. The contribution due to the large rj' mass dies off too rapidly, 
as seen from the asymptotic behavior. In this sense, it is very hard to measure the t]' 
characteristics from an x-space analysis of the correlation functions. The x-space analysis 
of the topological susceptibilities offers a better probe [^, although on the lattice there 
may be subtleties related to the definition of gluonic sources. Finally, we note that in the 
presence of the connected parts of the correlator, no asymptote sets in within 5 fm due 
again to the spurious oscillations discussed above. 

We note that since x^, = n ~ Nc, the instanton-induced shift in the 77' mass 2x*Nf j P ^ 



N^, at variance with Witten's argument [^. This is not totally surprising, if we recall 
that the original instanton gauge-configuration A ~ 1/(7 ~ V^iVc. Also for n ~ Nc, we 
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have a fixed compressibility a^, ~ N^. However, when the density n grows, the instanton 
and antiinstanton system is no longer dilute, and one would a priori expect a phase change 
23|1 , whence a breakdown of the conventional large Nc arguments. The academic case of 
n ~ yields zero compressibility, with the quantum fluctuations dwarfing the instanton 
effects. 

5.3 Scalars 

In Fig. 15 we plot the (normalized) connected C'^{x) and unconnected C^{x) parts of 
the correlator in short and long dashed lines respectively. The solid line represents the 
sum of these two. As seen from Fig. 3d the contribution from the first three diagrams is 
non vanishing in the scalar channel. If we were to repeat the calculation leading to the 
unconnected part of the correlator in x-space we would obtain the additional term 

N^jli S{k,m) (^^c^Tr S{k,m)-2 J ^ ^^ Tr C\k,m) m)/ + m))^ (59) 

If we recall the definition for the unconnected correlator we expect this additional term to 
be amenable to the square of the condensate < ipip >. The decay constant and the sigma 
meson mass follow from the last diagram of Fig. 3d. They can be evaluated through the 
use of a similar expansion of A_ to A+. Specifically 



d'^k Ml 



(27r)4 4P2 



and the scalar mass M_|_ is given by 



f\Ml . 2,„ i^Mr.) - (1 - ^) + / (61) 



+ \2N^ ' ' J (27r)4 V )^ V ' J 8tt^ 

where V = k"^ + M|. Numerically, we obtain m„=640 MeV and /=109 MeV. The scalar 
mass is about twice the constituent mass of 2 x 310 = 620 MeV. This is generic of 
all bosonized interactions at the mean-field level {e.g. Nambu-Jona-Lasinio model). The 
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nearness of the quark-antiquark threshold is expected to yield a large width for the scalar- 
isoscalar. 

5.4 Vectors 

• p, K* and </> 

To leading order in the instanton density, the vector correlation functions for both the 
p and (f) do not acquire any unconnected part. (The vector correlation functions are just 
the correlation functions of two screened quarks.) This is expected, since the instanton- 
antiinstanton interaction acts primarily in the spin-isospin zero channel. Figs. 16 and 
17 show respectively, the behaviour of the p- , and 0-correlators versus x up to 2 fm, 
for a light quark mass of 5 MeV and a strange quark mass of 140 MeV. The squares in 
these figures correspond to the instanton simulations, while the filled circles in Fig. 15 
correspond to the cooled and quenched lattice simulations. 

The failure to produce correlations in the vector channel, while obvious in the p-space 
analysis, is implausible from the x-space analysis. In general, simple spectral guesses as 
used in the instanton simulations or lattice calculations for an x-space analysis within 1 
to 2 fms may be misleading. They cannot differentiate between cuts and poles within 1.5 
fm. At these distances it is difficult to reliably differentiate between poles and cuts (the 
pion-channel being an exception). A resolution of the two requires a careful analysis of 
the preexponents and the asymptotics, as we have discussed. 

In the case of the K* , it is clear that a contribution due to the mixing between the 
up (down) and strange sectors occurs in the unconnected part of the correlator. This is 
evident from appendix E, where we see that the coupling of the p and to the quarks 
vanishes identically, whereas in the case of the K* a contribution in 0{ms — my) arises. 
The possibihty of the unconnected part of the ir*-correlator being contaminated by the 
excitations of its scalar partner (in flavour space) Trf- is allowed. Having said this we 
display in Fig. 18 the behaviour of the X*-correlator versus x up to 2 fm, for a light 
quark mass of both 5 MeV (upper curve) and 10 MeV (lower curve), and a strange quark 
mass of 140 MeV. Again, the squares in this figure correspond to the instanton simulation. 

• Ai and Ki. 
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Similar conclusions apply to the axial-vector correlators, although the latter are con- 
taminated by pion and kaon excitations through their longitudinal parts. Generically, the 
nonstrange axial- vector correlator can be decomposed along the tranverse and longitudinal 
directions that consist of the Ai and vr, respectively: 

C,Ap) = (5,, - P,P.)C^{p) + P.PuC^ip) (62) 

From the p-space analysis, each contribution is well separated. contains solely a cut, 
while displays only a pole. Similar remarks apply to the strange axial-vector correlator 
Ki. Figs. 19 and 20 show the behaviour of the combination "iC^ + versus x in the 
Ai and Ki channel, respectively. The squares refer to the results of simulations using 
instantons. 

Since the longitudinal pole reflects on the pion pole, consistency with the pseudoscalar 
correlators requires that the pion properties (mass and decay constant) should be the 
same. The explicit form of the longitudinal part of the axial correlator reads 



i?+ (p) 

c;:Ap) = p,p^^N^ ' (63) 

where at zero momentum 

which is just f'^/2Nc. Numerically, we obtain from the axial-correlator /,r = 76 MeV, 
which is about 10 % off from the value of = 88 MeV extracted from the pseudoscalar- 
correlator. This point illustrates some of the systematic uncertainties introduced by the 
use of undistorted scattering states for the nonzero mode states around a single instanton 



or antiinstanton 13 



6 Baryon Correlators 

In the large Nc limit, a baryon is made out of quarks, and is believed to be a soliton 



23]. In our case, we will think of a nucleon as made of A^c = 3 >> 1 quarks. To leading 
order in l/N^., the nucleon is just three free streaming constituent quarks. In contrast to 
the meson case, the induced instanton (or gluon interaction) interaction between diquarks 
is subleading in l/Nf.. We note that the soliton case in this model was considered in |2^ 
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Generically, the baryonic correlators will be defined to be 



Rix)=z {T*Jb{x) JsiO)) (65) 
where we use for the nucleon and delta currents 

J^{x) = eabc {u''{x)C^y{x)) u'{x) (66) 

respectively. Using Wick's theorem, we can reduce the nucleon and delta correlators into 
(Minkowski) 



R'^ix) = 2e„,eea'fe'c'7'^75^^'='(x)7"75Tr {j,S''' {x)j,S^''' {-x) 
R^{x) = 3e„,,e,,,,e'5^'='(x)Tr {^^S''' {x)Y 8'^'='' (-x)) (67) 



In the free case, (|67D reduces to 



respectively. 

Fig. 21 shows the behaviour of three constituent quarks versus x. The two solid lines 
are for 5 and 10 MeV, respectively, the open circles are the results of instanton simulations 
and the full circles are those of quenched and cooled lattice simulations. Clearly, both 
simulations show attraction in the nucleon channel, which is very likely due to the fact 
that in the instanton model, the instanton induced interaction in a spin-zero isospin-zero 
diquark configuration (qq)/z'o is attractive. This follows from the large attraction already 
observed in the spin-zero isospin-zero quark- ant i- quark configuration {qq)j^Q by crossing. 
This attraction is, however, an order of magnitude smaller than the attraction in the pion 
channel. Whether these interactions can result in a pole remains an open question and 
requires a more detailed analysis. Indeed, the screened quarks amount to a mass of about 
940 MeV, which is close to the empirical value of the nucleon mass. 

Fig. 22 shows the results in the delta channel. From this, we conclude that the 
present simulations cannot distinguish between a cut and a pole in this channel. In fact, 
it is very unlikely that a dilute instanton gas can yield binding in decuplet channels, since 
the instanton induced interaction is usually non-existent in these channels. 
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7 Gluonic Correlators 



The present construction allows for a convenient analysis of correlation functions involving 



F ■ F and F ■ F reflecting on the scalar and pseudoscalar glueballs in the model ||2T|, [17 
In the quenched approximation, these correlators are ultralocal and given by our choice 
of the measure (^. Through the identification 

1 



327r2 

the scalar gluon correlator reads 

327r2 ^ ^ 327r 



F ■ = (n+ + (69) 



Cff{x -y)= < T*^F\x) ^F\y) >conn.,N,=0 



al5\x-y) (70) 



Also, through the identification 

1 



F ■ F{x) = {n+ ~ n-){x) (71) 



327r2 

the pseudoscalar gluon correlator reads 

Cff{^ -y)= < T*^,FF{x) ^,FF{y) >^,=o 

= < T*(n+ - n'){x) (n+ - n~){y) >Nf=o 

= X* S\x - y) (72) 

In (|^),the glueballs in the quenched approximation carry infinite mass and zero size. They 
act as heavy sources. 

In the presence of quarks, the glueball sources mix. The mixing is of order In 
the scalar channel, 

Cff{x — y) = <T* [^n^ + rr){x) — {{n^ + n^){y) — > 

= al(^ 6\x -y) + 2Nfal < T*7io{x)My) > ) (73) 

where the unconnected correlator in the Aq scalar channel is 



< no(x)no(y) >= 4r fdk e^'^^""*') ( —-^ r + , r ) (74) 
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The large separation behaviour of the above result follows from section 5 with the pole 
mo = 640 MeV as the mass of the scalar-isoscalar. Because of the mixing, the fall-off is 
dictated by the scalar-isoscalar masses. Fig. 23 shows the plot of the scalar correlator 
Cff{x) (minus the ultralocal term). From (|73| ) and ([f4D, the compressibility takes the 



form 



a^ = l.^<(^J (fz{n+ + n_ - >^ y (75) 



where h is given by 



IIN^ 2Nf , , 

b = - —^a+ (76) 



with 



Numerically we find a+ = 1.22, which is to be compared with a;+ = 1 in the QCD trace 
anomaly. This is only suggestive, however, since the two calculations are totally different 
in spirit. Ours is classical, while in QCD it is quantum. 
In a similar way, we have in the pseudoscalar channel 

Cppix -y) = < T*{n^ - n'){x) (n+ - n'){y) > 

= X* ( ~y)~ 2NfX. < T*r]o{x)7]o{y) > ) (78) 

where the unconnected correlator in the Aq pseudoscalar channel is 



The rotation angle 6{k) along with X± are defined in ( pOD and (|^), respectively. The 
fluctuations in the pseudoscalar gluonic source fall off with a rate that is given by the 
lightest mass (the rj in our case) = 557 MeV. Fig. 26b shows the plot of the scalar 
correlator Cpp{x) (minus the ultralocal term). Let us now evaluate Cpp{x — y) using the 
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pseudoscalar field decomposition, in which the quadratic part of the action is diagonal. 
We obtain 



r /I 1 ^/ 1 \ -1 

Cj,p{x-y)= dk e'''^~y^ — + TTTT E IT^, ^ (80) 



One should not be alarmed by the two different expressions for the pseudoscalar glu- 
onic correlators (|78D and (|80|). Using the two relations (denoting A = A_(A;, m,m) and 
As = A_(fc,ms,ms) ) 



A+(A:)A_(A;) = 4AA, + ^^(A + 2A,) (81) 



and 



A+(A;) sin^ d{k) + A_(A;) cos^ d{k) = -(A + 2A,) (82) 

3 



we can easily rewrite (|80D to (iTq). The mixing causes the topological susceptibility to 
decrease. From (|80D we have 

X = ^^<{ fd'z{n^-n-){z)y >= ^ (83) 



and vanishes for any quark mass going to zero. The topological charge is totally screened 
in the chiral limit. 

8 Nucleon form factor 

All hadrons are characterized by various form factors, each of which carry information on 
the various charge and current distributions. In this part, we show how various nucleon 
form factors can be analyzed in 1/Nc, thinking of A^c = 3 >> 1. In this section, we will 
distinguish between purely gluonic form factors G(x) ~ F^(x) , FF{x), a^^uF^^F^ , ... and 
fermionic form factors F(x) = ip'^ Aip, where A = 7 T is a spin-flavour matrix. Mixed 
form factors M(x) = ip'^a^yF^^ip, ... can be obtained in a similar way, although they will 
not be discussed here. 
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• Gluonic Form Factor of a Constituent Quark 
Since the model lacks confinement, the nucleon form factor receives contribution from the 
unconfined constituent quark states. This is represented in Figs. 26(a) and (b) as 
insertions. To leading order in l/Nc-, these contributions are either direct as shown in 
Fig. 26(a), or meson mediated (Figs. 26(c) and (d)). The constituent quark gluonic form 
factor is defined (when x — oo) as 

^G{k') (T>(|)V^t(_|)) = (T>(|)G(fc)^t(_|))_ (84) 

with 

G{k) = I dye'^-yGiy) (85) 

Throughout, we will think of Fq as matrix valued (here in spin space), so that various 
components of the form factor can be extracted by proper tracing. To leading order 
in 1/Nc, the form can be readily evaluated using the bosonization results developed in 
appendices C and D. The result is (x oo) 

FpFik^) S{x,m) =i Jdpe'P-''^Mp_Mp^C{p^) 32n^ {Ti,{-k)a{k))C{p+) (86) 

for a scalar gluon insertion, and 

Fppik'') S{x,m) =1 J dpe'P-^^Mp_Mp^C{p.) «327rS(vrp.(-A;)x(A;))C(p+) (87) 

for a pseudoscalar gluon insertion. Here p± = p ± k/2. The mixed spin-gluon matrix 
element can be discussed using similar arguments. The expectations in ( PB[ - P7D involve a 
Gaussian integral over the effective bosonic fields, with quadratic actions as discussed in 
Appendix C. After integration, the results are (x^ oo) 



FpFik'') S{x,m)= + t I dpe'P-''^MpMp+Cip^)Cip+) 



327r2 / n 1 ^ 1 \-i 



2NcA^{k, m, m) xn^al 2Nc A+(/c, m/, ruf) 



and 



Fp^iP) S{x,m)= + t J dpe"P-^^Mp_Mp+C{p^)^^C{p+) 



327r2 / 1 1 



X 



- + T77rE . J (89) 



2NcA^{k,m,m) \x* 2Nc^ A_{k,mf,mf) 
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From our numerical analysis of section 3, the constituent quark propagator S{x,m) 
shows a rough scaling in the window < x < 2.5 fm, with Mo ~ 300 — 400 MeV, but 
then oscillates for x > 2.5 fm, due to non-analyticities. In the window — 2.5 fm, 



It would be interesting to see how the present form factors ([88| - |89D with (|90|) compare 



with simulations in the range < x < 2.5 fm. This is only indicative, since the channel 
is contaminated by spurious oscillations for x > 2.5 fm. 

The large x separation provides for a way to select the constituent quark on its "mass- 
shel"0, hence the analogy with the Minkowski definition of the form factor. We can also 
define a totally "off-shell" form factor by considering (8^89) for finite x and integrating 



X over V4. In this way, one obtains off-mass shell form factors with zero-momentum 
constituent quarks. For k = 0, the results are 

Fff(O) = 2iV,A_(0,m,m) 2iVe ^ A+(0, m^, m^) ) ^^^^ 

and 

" 2iV,A_(0,m,m) ^ 2iVe ^ A_(0, m;,m;)) ^^^^ 
To leading order in l/N^, the "off-shell" scalar form factor reduces to 

F^f(O) = (93) 



where a = 2a+ ~ 2.45. We note that ( P3| ) differs by almost a factor of 2 from its "on- 
shell" analogue with a = 1, as argued from a QCD low-energy theorem based on the trace 
anomaly |26|, For the pseudoscalar form factor (|9^) the result is F^^(0)=0.44, which 
is to be compared with the gluonic part of the "on-shell" value of the axial-singlet form 
factor g^, as determined from the U(l) anomaly (^ in the constituent quark state 

< r*^t(|) ^ J + jfjdz Tr; m^t^^^(^)^^(_|) (94) 



^This is, of course, suggestive in Euclidean space. 
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The mass term in (|9^) involves the Ua(1) form factor in the constituent quark state. In the 
last few years, efforts have been made to understand the data from the European Muon 
Collaboration (EMC) |^6|. One of its remarkable results has been to yield a small 

value for the singlet axial coupling constant g% = 0.13 ± 0.24. The result obtained above 
"off-mass shell" seems to be close to this value. The approach described here, provides 
some insights from a instanton vacuum model, to the effective approach discussed by 



many p8| , ^ In fact, the modified bosonisation scheme discussed in Appendix E, is 
very close in spirit to these models. A comprehensive discussion of all these issues goes 
beyond the scope of this work. 

• Fermionic Form Factor of a Constituent Quark 



The fermionic form factors can be analyzed in the same way as the gluonic form factors. 
The mechanism consisting of insertion is shown in Fig. 27. In Fig. 28, we show the 
leading contributions to the mesonic form factor to order l/Nc- Fig. 28(a) counts the bare 
charge, while Fig. 28(b) involves a typical meson-exchange with non-local form factors. 
Generically (x — > oo), 

FA{k') (T*^(^)V^^(-|)) = T^i^) ^^\m ^H~))conn (95) 

Parametrizing all meson fields by = 'j'^n^, where 7^ = (1,75) (S> T, yields to leading 
order in 1/Nc 

Fa(A;2) S(x,mj) = - J dpdqe'^P^^/^^'' ^pC{p, rUfhCiq + k, nig) ^ M^+u A 

X V^C(g, mg) 7 C{k -p,mf) ^fM^^^-K^^ip - q - k) n'^ {q + k - p)) 
+ f rfpe^(P-^/2).. ^^c{p,mf)jC{p - k,mg) ,/M^k{n^\k) n^^{-k)) 



X Ti[A^ M,C{q,mg)^C{q + k,mf)^M,+k) (96) 

where summation over flavour g is understood and / is the flavour of the quark being 
probed. 7 is 1 (75) for the scalar (pseudoscalar) sector. The first and second terms of 
(p^) are displayed in Fig. 28(a) and 28(b), respectively. The expectation value involves a 
Gaussian integration over the measure derived in Appendix C and can be evaluated for 
arbitrary momentum q. In the scalar sector 

1 1 



^ ^ ^ ^ 2N,AJq,mf^mf) \ 



2NcA+{q,mf,mf) V 2A'cA+(g, m/, m/) 
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^ n,a^,^^2N,A+{q,mg,mg)) ^^^^ 
For the pseudoscalar case, we replace A+ by A_ and n/{n^al) by 1/x- The case where 
1 = IfjL (vector form factors) and 7 = cr^j, (tensor form factors) can be analyzed similarly. 

It is interesting to note at this stage that most of these form factors may be used to 
assess the strength of the meson-constituent quark interaction in some constituent quark 



models, as recently discussed by Glozman and Riska ||9|. When couched in the 1/Nc 
framework, the present analysis provides some rationale for their successful phenomenol- 
ogy- 

As in the gluonic case, we can investigate the "off-shell" limit of the form factor at 
/c = 0. Using (|96|) for fixed x, integrating numerator and denominator over the entire V4, 
and taking the k = Q limit, yields 

FA(0)^(0,m)= - M^C{0,m)-f Jdq Mg 

X (c{q, m) A C{q, m)7 - Tr (A C(g, m)^ C{q, m)) ) C{0, m) (98) 



where all momenta are taken to be zero. Numerically, the meson-meson expectation value 
(7rs(0) 7rs(0)) in the scalar sector is 0.69 fm ~^ and 0.92 fm for the up (down) and strange 
quark, respectively. The same applies for the pseudoscalar sector, where (vrps(O) TTps{0)) is 
4.68 fm and 6.69 fm for the up (down) and strange quark, respectively. In short, 
formula (|98D along with the numerical values for the meson-meson expectation value, will 
serve us as a check point when numerically generating the values of Fa(A;^) using (pBI) . 



• Form Factors from loffe's Currents 
If we were to think about the nucleon as made out of three constituent quarks, then the 
nucleon form factor follows from the additive constituent quark picture. When simula- 
tions are performed, however, it is customary to use loffe's currents (^) for the nucleon. 
This results in some non-trivial combinatorics and folding of the single constituent quark 
propagators, as we now explain. Let J^{x) be loffe's current (|66D . Then, the nucleon 
form factor reads {x —>■ 00) 

F^iP) (rV^(|)J^(-^)) = (TV^(|)0(fc)jJ(-|)) (99) 

where O = G, F, which are short for the gluonic and mesonic insertions discussed above. 
Typical diagrams for mesonic insertions are displayed in Fig. 29. The term (L.H.S.) 
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multiplying F^^k"^) in the left hand side of ( p9D can be readily reduced to give (|67D. The 
right hand side (R.H.S.) takes the form 

R.H.S. = +6 (7^75(95(0;; fc;m)7p75)^^Trs (7^S'(x,m)7pS'(-a;,m)) 
+6 (7^75S'(a;,m)7p75)^^Tr, (7^05(x; k;m)-fpS{-x,m)) 
+6 (7/.755' (x, m)7p75)^^ Tr^ (7^5* (x, m)-fpOS{-x, k, m)) (100) 

where OS{x] k; m) follows from the right-hand side of (p6|) and (pT]) for the gluonic inser- 
tions, and, (^) for the mesonic insertions. 

It would be interesting to see how ( |100D compares to actual simulations. As noted 



above, the actual constituent quark propagator oscillates at distances larger than 2.5 fm. 
Hence, a true asymptotic form factor may not be reached in this model for the nucleon. 
In the region < x < 2.5 fm, the constituent quark propagator seems to be damped 
following the behaviour described in (|90|) . Using this behaviour, the left hand side term 



in (|99|) reduces to 

L.H.S.= fl^,/;^e-«»Y 



xQ [ipi,OS{^ + l)7p75j^^ Tr, [^p{^ + lhp{~^ + l)j (101) 



while the right hand side reduces to 

R.H.S.= 6 f^,/^e-H' 



X (^(7M75C5(x;/c;m)7p75)^^Tr, (7^(^ + l)7p(-^ + 1)^ 
+ (7^75 (^ + l)7p75)^^ Tr, (^7pO>S(x; fc; m)7p(-^ + 1)) 
+ (7;.75(^ + l)7p75)^^Tr, {^p{jt + l^^OSi-x; k;m)) ) (102) 
Numerical results for the resulting form factors will be given elsewhere. 

9 Discussion 

We have analysed the mesonic correlators in a random instanton gas in momentum space 
using bosonization techniques, and, in coordinate space by performing direct Fourier 
transforms. Our starting point was a grand-canonical ensemble of instantons and antiin- 
stantons, where the t' Hooft vertices play the role of "fugacities" . The momentum space 

31 



results are in agreement with the original analysis in both the massless and massive 
cases 0]. Following 't Hooft's suggestion, the resolution of the Vj problem follows by 



assuming that the topological charge is screened [T^, with a finite screening length 
(non-zero topological susceptibity) . This effect is leading in IjNc counting and results in 
a contribution of order A'"^ to the r/' mass. Without this effect, the r( would be degenerate 
with the 77. 

We remark that a non-vanishing topological susceptibility should not be taken for 
granted 0. In the present case, it follows directly from the use of instantons and antiin- 
stantons in a singular gauge. A check would be to repeat the analysis using instantons and 
antiinstantons in a regular (non-singular) gauge, or, carry out cooled lattice simulations 
with free boundary conditions. 

Our x-space translation of the p-space correlators shows that the results of simula- 
tions using either a large sample of instantons and antiinstantons in four dimensions, 
or quenched and cooled lattice gauge configurations, are in agreement with the Fourier 
transformed analytical calculations within the reported range of (0-1.5) fm. The recent 
analysis carried out in Ref. |^ for two flavours differs from the bosonized results [|I^ [], 
hence our analysis. 

We have shown that the running quark mass causes the quark propagator to oscillate 
at large x. The oscillations are larger for larger quark masses and affect most of the 
correlation functions at large distances. These effects are spurious and reflect on the lack 
of confinement in the model. They are easily subtractable in a p-space analysis. They 
are harder to track down in an x-space formulation. The extent to which these spurious 
modes impact on the subtracted results is presently unclear. 

We have shown that, while the asymptotics of suitably subtracted correlators yield 
pseudoscalar masses that are accurate to within a few per cent, the non- asymptotic read- 
ings could be as inaccurate as 100 %. From our calculations, the subtracted and rescaled 
correlators show good asymptotics between 2 and 3 fm. The non-rescaled correlators do 
not show any reasonable asymptotics even up to 10 fm. This point merits further scrutiny 
in lattice calculations. 

The bosonized results show that while it is possible to infer the existence of light 



^Eq. (58) in Ref. |p5| relies on a resummation of the quark propagator Eq.(57) which is vahd only for 
zero quark mass. 
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pseudoscalars in a dilute instanton gas, they do not seem to support the appearance of 
bound vectors. We have exphcitly shown that the resuhs of simulations are consistent 
with the presence of just screened quarks in these channels. We have noted that the use of 
schematic poles and cuts to analyze the x-space correlators in these channels would have 
imphed otherwise. Due to mixing between the octet and singlet pseudoscalars, we have 
found it difficult to extract the t] and t]' masses from the x-space analysis. The extraction 
is straightforward in the p-space analysis. 

We have presented a simple analysis of the baryonic correlators in both the nucleon and 
the delta channels. The attraction seen in the nucleon channel is expected from general 
arguments. In this channel, however, it appears to be difficult to identify a nucleon mass 
without going to the asymptotics, since three screened quarks already yield a mass of 
the order of 940 MeV. This may cause the nucleon to unbind, although soliton-inspired 
calculations with constituent quarks seem to suggest otherwise [^. In this respect, it 
would be interesting to repeat our analysis by including diquark fields. The results of 
simulations in the delta channel are also consistent with three constituent quarks. A 
dilute instanton gas does not induce correlations in the decuplet channels. 

Using loffe's current for the nucleon, we have worked out various gluonic and mesonic 
form factors "on- and off-mass" shell, to leading order in l/Nc- The form factors are 
sensitive to the the three constituent quark cut. Moreover, the appearance of spurious 
oscillations in the single constituent quark propagators causes the form factors to be ill- 
defined for point-to-point separations that are larger than 2.5 fm. In the region < 
X < 2.5 fm, some estimations have been made that would be of some interest for future 
simulations. The analysis of the nucleon form factor presented in this work could also 
be extended to other mesonic and baryonic channels. It also provides insights in to some 
recently used constituent quark models [^]. 

The fiuctuations in the number sum and difference of the instantons and antiinstan- 
tons relate directly to the scalar and pseudoscalar glueball correlation functions. In the 
quenched approximation, the glueballs are infinitely heavy and stable. In the unquenched 
approximation, they mix with their scalar and pseudoscalar counterparts and decay. The 
mixing and decay are of order l/N^. 

The overall agreement between the instanton simulations and the present analysis 
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within 1.5 fm shows that a random set of instantons and antiinstantons that is suitably 
stabihzed in the infrared is well described by gaussian fluctuations over a mean field 
solution. The mean field solution follows from a simple bosonisation scheme. It also 
shows that constituent quark models with dynamically generated masses, e.g. Nambu- 
Jona-Lasinio model, are also likely to give similar results provided that chiral symmetry 
is dynamically broken. In all these models, however, the subtle issue is that of confinement 
with its impact on large distance asymptotics and form factors. 
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10 Appendices 



10.1 Appendix A : Generating Functional 

In this Appendix, we provide the necessary details for the derivation of the generating 
functional ( |15|) discussed in the text. Although these calculations were extensively used 
in establishing the results of Refs. |jl9[ , they were never published. We start by evaluating 



the color averages occuring in the 't Hooft determinants (|T]) for Nf = 1. For convenience, 
we will use the shorthand notation d^k/ (27r)^ — > dk and d'^x —>■ dx when integrating out. 
If we denote by 

9^{z) = { J dxij^S^'(P^{x -z) j dycp^^y - z)S^^,lj{y))up (103) 
then its Fourier transform reads 

e^{z) = jdkdi e"'(^"')" e^{k, I) (104) 

with 



Averaging over the color group, we obtain 



e-(kj) = ((i - ^) of (i - 9^));fe«W (106) 

where (j)'{k) is the Fourier transform of the fermion zero mode profile, and is given by 



(t>\k) = 7Tp^-^{Io{z)Ko{z) - h{z)Ki{z)),=,p/2 (107) 
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With the use of ([T^) , the partition function (0) takes the form 



X exp ^/rf.log(l-^^^(.)) 



X 



exp i Jdkdl P^{k,l) (n^{k,l) - e^{k,l)) (108) 



where the integral in the last exponent is performed in both variables k and / of the bifocal 
auxiliary fields. The field vr^ is eliminated using the mean field equation 



ri p 1 p i(^k l')z 

-'^^(^■O-j/''^ r^^E^^ (109) 



2im 



For Nf > 1, the auxiliary fields vr^ and are Nf x Nj valued along with the average 
6^ ( 103 ) entering the 't Hooft determinants (|l]). As a result, additional traces over flavor 



indices will be needed. With this in mind, the previous results can be generalized in 
a straightforward way. For m = diag(mi, mjy^), the result is (after absorbing in the 
measure a term in the size p to have a dimensionless argument in the log) 



where S ^[P^,P ] is given in the text (p^. At the saddle point P^ = P, and S is the 
quark propagator in the external background P such that in momentum space 

S{k,l) =6\k-l) S{k,m) (111) 

with S{k,m) written down in (^). From (|110|) , the partition function (|15]) follows after 
integration over the fermionic flelds. 
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10.2 Appendix B : Quark Condensate 

In this Appendix, we will show that (^) follows from an exact derivation using the 
standard definition prior to the bosonization procedure. Following the method used in 



Ref. [19 1, the partition function (0) can easily be written as (ignoring fluctuations in the 



density and switching off the sources) 

Z={j Vxjj^ViPe-S^^^''^) (112) 

where in the one flavour case 

= S,' + ^So%4So' (113) 
The (Euclidian) quark condensate follows as 

<^t^>= J—{fv^l:^Vij [ ^/^Ve"/^^^"''^) (114) 
where averaging over all pseudoparticles is understood. Specifically, 



,..t..,- {^rS{0,m) detj-S-^)) 

{det{-S-^)) ^^^^^ 



Introducing a set of Grassman variables for the pseudoparticle ensemble, the partition 
function (|112|) reads (sum over J, J understood) 



Here, T is the kinetic part of the overlap matrix [|13| [l9| and the integration is over 
fermionic fields if), ip^ and Grassman variables xii xl? where / is an integer that runs over 
all the instantons and antiinstantons in the ensemble. Similarly, the condensate 



< 



From the two formulas above, it is easy to show that 



>= <^t^ >= --1^1^ (118) 
Va dm 
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We now shift the fermion fields according to ip ^ ip + i(j)iXi and ^ + ^xl^l (sum 
over / understood). We then expand x ^iid around the respective classical solution of 
the shifted action. The remaining integral ( |117|) has now a Gaussian form in x^X ^^id can 



be performed. We obtain 

< ^ V > = 77^ ( / Vij^Vi; (-2«m)^ e" / '^'^"''^ 

1/4 Z/ J 

y lyimy j im J 

Rewriting the pseudoparticle sum in the exponent as a product over /, and noting that 
only the first two terms in the Taylor expansion contribute, we can easily perform the 
color group average to yield 



where 0^{z) is given in appendix A. As in |T^, we assume a sufficient amount of coarse 



graining so as to rewrite the product over / with the result 



Va Zj j 



+ ^ jdz^-^^}$l^\ et^^Mi-^e±W) (120) 
^ ^ 2imJ 1 - e^{z)/2im / ^ ' 



where the dot on 9^ {z) indicates the derivative with respect to m. The functional integral 
above can be evaluated exploiting the same bosonisation scheme used for the partition 
function Z 

<^V> = -i- /p^tp^I)P±(-2zm)^e-/'^'^''[^+'^"l'^ 
Vi^Z J 

X ^{-^+^-^Sd^PHz)) ^t^ + 1^^ _ 2*) P^{z)^ (121) 
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where S ^[P^, P ] is given in the main text (|T^. At the saddle point 



<^t^> = J_ fv^^V^(^J!:^fe- I^'S-^^^-Ni-N+4mP) 



X 



J (-1 + 2SoS-^) ip - AiVP^ (122) 



where (aside from rescahng Mk{m)) the quark propagator 5* is written down in momen- 
tum space (0). After performing the integral and properly rescaling P, we recover the 
expression (|3TD quoted in the main text for the condensate (in the chiral limit). This 
result is expected, since to leading order in l/Nc, the determinants in (|115|) cancel out, 
after factorization (quenched approximation). 

10.3 Appendix C : Gaussian Approximation 

In what follows, we give details leading to the Gaussian approximation in the partition 
function. We can repeat the steps performed in appendix A with the constraint = 
= n/2 now relaxed and the parametrization 

P^{z) = P + n^{z) (123) 

A few comments are in order. In the equations above, a{z) and xi^) respectively represent 
the scalar and pseudoscalar glueball sources. The field vf contains pseudoscalar and scalar 
excitations and will be discussed further below. Following appendix A, the auxilliary field 
TT^ is eliminated using the mean field equation 

-tP^(k,l) = — dz- -f— (124) 

^ ' 2 J 2im 1 - 7i^{z)/2im ^ ' 

Along with the contribution from the measure /i(?T-+, n_) (|^) the bosonized effective action 

reads 



S^Q= - N,Tt log S-^ [P+, P-] - 2 JdzTTfm(^P+{z) + p-{; 

+ Jdz n±{z) Tr^log (^ ^^^^ ^ + Jdzn±{z) logNfl + So (125) 
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where 

Sg= + ^ I dz {n+{z) - n_{z)Y 

+ —Jdz (n+(z)+n_(^)) (^log +^ ^ ^ ^^^^^ 

The trace (Tr^) is in flavour space and the trace (Tr) is over flavour and Dirac indices 
with an integration over momentum. 

• Gluonic contribution 

Let us first turn our attention to the last three terms of the effective action 



Sg 



+ Jdzn±{z)TTflog\-—^ 
+ Jdz n±{z) log Nf\ 



+ 



2x* 



dz {n+{z) — n_{z)y 



n+iz) + n^iz) 



n 



(127) 



Using the saddle approximation in the scalar glueball source a{z) fluctuations, we obtain 

=nexp ^^logNflW — ) (128) 

As flrst discussed in and later in [0, the distribution of the fluctuations in the num- 
ber densities n±{z) is Gaussian (exact) in xiz) with a width given by (^. The distribution 
is logarithmic in the sum oi^z) and Gaussian (approximate) in the large Nc limit with a 
dispersion relation given by (|^). Along with the saddle point decomposition of the bifocal 
auxiUiary fleld = Pe^'''^'^'^{l + Trje^^*'''"'/^ , we obtain 



Sg [P^,ni] = [7r,,p,] + Sg^ [vr,,^,, a, x] (129) 



^This parametrization is reminiscent of the action being invariant (for massless quarks) under global 
axial transformation with the subscript s and ps respectively standing for the scalar and pseudoscalar 
mesonic excitations. 
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Adopting the nonet decompositions Tips = Xqtjq + J2 -^aT^ps and vr^ = XoTigfi + Xair'^, the 
term Sq^ ['^s,ps] contains mesonic fluctuations only and reads 

4^ [^s,ps] = -y /rf^ Tr; (n'.iz) - nliz)^ + jdzN^ {x*vl{z) - a>Jo(^)) (130) 

We point out that the term Sq^ ['^s,ps\ should be put in concert with the first two terms 
of Ses[P^ ,n±\ in order to obtain the total mesonic contribution to (|125|) . 



/ox 

The last term 5*^. [T^s^ps, cr, x] involves mixing on the one hand between the isosinglet 
scalar and the scalar glueballs, and, between the isosinglet pseudoscalar and the pseu- 
doscalar glueballs, on the other hand. 

S'§^ [Ks,ps,(r,x] = + Jdz ^(^xiz) +ix*^J^fVo{z)^ 

+ jdz^{a{z)+al^fTis,o{z)^^ (131) 

• Mesonic contribution 

Performing a Taylor expansion of around the saddle point P in the first two terms 
of S'eff[P^,n±] (|125D along with [ vr, „ J the total mesonic contribution reads PI 

+ N,ldk7rls,{k)A^{k,mf,mg)nf^{-k) 

+ Jdk {vo{k)^Vo{-k) - vr.o(A:)y^vr,o(-A;)) (132) 

where the saddle point approximation leads to an integral (gap) equation in P{mf) for 
each flavour / 

4iVc /■ „ .X X 2P(mf) 



n 



JdkA{k;MkP{mf);mf) = 1 - 2m/ (133) 



^The sum over flavour indices / and g is understood. 
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After rescaling the constituent mass according to M^P —>■ with P{mf) = raA(mj)/2, 
we obtain the gap equation (^OD in the text. We define below the various quantities 
introduced in the mesonic action ( |132| ). In momentum space, we write the inverse quark 
propagator in the background of instantons and antiinstantons as 

{k\S-\P)\l) = 6{k-l)S'^{k,m) (134) 

and 

S~' {k, m) = - .m) (^^ - z - j (135) 

The coefficient A [k; M^; m) appearing in the gap equation is given by 



A [k; Mk, m) = — ^ (136) 



Except for the isosinglet scalar and pseudoscalar, the inverse meson propagator in the 
background of instantons and antiinstantons, apart from the factor f'^/ANc, can be iden- 
tified with A-|-(A;) appearing in the quadratic part of ( p.32| ) and reads 

A±(fc,mi,m2) = ^-2 j dqiA^A^T {qi-q2) BiB^) (137) 

where we have set gi,2 = 0' ± Mi =Mq^{mi), Ai = A (gi; Mi; mi), Bi = B {qi, Mi, mi), 
mi being one of the quark masses in SU (3) flavour space and B is given by 



t.2 1 

B (k; Mk] m) = — ^ (138) 

^ ' ^ Mkk^ + {m-kyMkf 



In what follows, we will always consider the rescaled constituent mass Mkimj 



• Bosonized partition function 
To be thorough, let us exhibit the bosonized partition function utilized in evaluating the 
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various (mesonic, baryonic and gluonic) correlators. To this end, 



Z = ypvr.^p, Dct e"^^«f^^'''^''''^' (139) 



where the bosonized action follows from regrouping terms in Sq (|129|) and S'meson ( |132|) . 



5'cfr[vrs,ps, x] = Sf^^, ^, 0] + S^s\t^sA + 'S'efr[o-, x] (140) 

with 

'S'iff [0,^,0] = -N^llxlogS-^[P) - ^ ^ NHx^Xog — - Wm^P{m^) (141) 



along with the mesonic part of the effective action 

5'efT[vr,,p,] = - N,Jdk7rf'^{k)A+{k,mf,mg)7if{~k) 
+ N,Jdknfj{k)A^{k,mj,mg)nf^{-k) 

+ Jdk (^r^o(fc)^r/o(-A:) - 7r,o(A:)^7r,o(-A:)) (142) 
and the gluonic part of the effective action 

+ jdz^^[a{z) + al^f7:,^,{z)^' (143) 



We are now in a position to evaluate the correlation functions of interest. 

• Connected meson correlator 
From the expression of C^{x) (^) in the text along with the bosonized partition function 
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( |139|) , we easily find that to leading order in Nc (the trace being over flavour as well as 
Dirac indices) 



C°(x) = —NcTt(^ S{x,m) 7 S{—x,m) 7 



(144) 



along with its p-space version C°(p) quoted in ( P7| ) of the main text. 

• Unconnected meson correlator 
From the expression of C}^{x) (^) in the text, we need to examine the term Tr '~fS{x, x, P^) 
in the integrand (the trace being in flavour, color and Dirac space). With the shorthand 
notation 



vr, 



ps 

2~ 



vr = 7r.,- — + ^75ivr,. + -vr.vr,. + -v.j 



(145) 



we can write 

Tr jS{x, X, P^) = jdkdl e*(^-^)^ Tr 75(A;, I, P^) (146) 
where the relevant term in the large Nc limit is given according to 

S{kJ,P^) = S{k,m) 6{k-l) + tJ^^C{k,m)n{k-l)C{l,m)JMi 



- ^/ATkCik^m) (^Jdqn{k-q)^q(^l 



(147) 



The coefficient C{k,m) is given further below. Defining 

C}^{p,q)=6{p-q) (C';(j9)+C^'(j9)) 

where 

S{p) 



C'ip) 



V-Ks^ps Va Vx exp ( - Scq[ks,ps, o", x]) 



(148) 



X Tr 7S Tr 7^ - 2Tr -^jMhC ^ {Bl^ + iA) ^ C J (149) 
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Here, the trace carries an integral over momentum. The meson field vr is left inside 
the trace so as to refiect its bilocal character in the momentum. It is clear that C}^{p) 
vanishes identically except in the isoscalar singlet vr^ o- From a diagram approach, it has 
two unconnected closed fermion loops with possibly the tTs q being emitted within one 
loop (Fig. 3c). Further discussion will be presented in the text regarding this term. 

We are therefore left with the second term contribution in C}^{p). This term amounts 
to a propagating meson from one closed fermion loop to the other (Fig. 3d) and reads 

X exp - Sc{i[7is,ps, cr, x]j (150) 



where 



R^{p,mf,mg) = Jdk ^ MrM2Ti{iC{ki,mf)n{p) C{k2,mg)^ (151) 



with ki^2 = k±p/2, Ml = Mk^ijUf) and 



Cik,m) = Sik,m)(l-'-^^ (152) 



If we redefine so as to extract the meson field, we have 



with 



Ry{p,mf,mg) = R'^{p,mf,mg) ± R'^{p,mf,mg) (153) 



(P, ^1,2) = / ^V^MiM^Tr (tCi^C^) (154) 



and 



h2 
'^1,2 



Ci,2 = T^il-Ai,2)+tBi,2 (155) 



Performing the functional integral, we obtain (pS]). As an example, we exhibit the case of 
the mixing singlet 7]q and octet r/g- The specific fiavour character of R^ follows from the 
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effective action ( |142| ) for the nonet decomposition. For instance, the pertinent terms R^^ 
( [15 1|) for the r^o and r/g correlators are given by 



dk^MiM2TTf (75Ao,8C(A;i, m)^5C{k2, m)/s:(p)) 



(156) 



Explicitly, for the unconnected parts of the tiq and r^g correlators, we obtain 



CLM = -f 1-^^0,8 Ro,8(-p)r/(-p)77(p)Ro,8(p)e^^/^[^]'' (157) 



where the partition function Z in the denominator contains only tjq and r^g. For the rjQ we 
have 



Ro(p) 



/ lR(p) + lR^{p) 
[ ^-^(R{p) - R,(p)) 



(158) 



and for the r^g 



R8(P) 



V |/2(p) + |i?,(p) 



(159) 



with R{p) = R-y^{p,m) and similarly Rs{p) = -R^g(p, m^) where m and are the up 
(down) and strange quark mass respectively (|151|) . 



(160) 
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10.4 Appendix D : Nonet Decomposition 

Taking the partition function (|139|) derived for Nf > 1 with the same decomposition for 
the meson fields vr^^ps, we have as the mesonic effective action ( p.42|) 

Ses[^s,ps] = - Nc Jdk 7c{^{k)A+{k,mf,mg)nf{-k) 
+ N.Jdk 7ifj(k)A^{k,mf,mg)7if^{-k) 

+ Jdk (^^,{k)^r^,{-k) - vr,o(fc)^vr,o(-fc)) (161) 

where / and g are flavour labels (f is not to be confused with the pseudoscalar decay con- 
stant). Using the decomposition Tips = J2k=o^k "^k, the kq and Kg excitations contribute 
to ^eff in the form 

Ses[vo,V8] = N,Jdk^{k)[A{k)]7]{-k) (162) 



where f/(p) = {vo{p) , Vsip)) = f (i^oip) , i^sip)) and 



[A{k)] 



I i(2A 
V 



A,) + 



Nc 



^(A-A, 



¥ (A - A,) \ 
f(A + 2A,) ) 



(163) 



with the shorthand notation A = A{k,m,m) A^ = A{k,ms,ms) with the strange quark 
mass inserted. The contribution x*^f/^c follows from the singlet mixing with the topo- 
logical fluctuations through the measure (^. At low energy, and to leading order in the 
current mass 



SAVoM = + jdk-fk\ril + ril)+ j -{2NFX*)ril 
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f 1 — / 4: 2 2 4 

- J dk-{i)ij) i^r]l{-m + -nis) + ril{-m + -m, 

f 1 — 4\/2 

- J dk-{^tp)r]oris—^{ms-m) (164) 

Note that 5'efr[?7o5 "/^s] ~ ^c- The above result yields the GOR relations for the singlet and 
the octet, if we were to drop As is well known, the GOR result is badly violated 

in the singlet channel by the axial U(l) anomaly. The latter is carried over by local 
fluctuations in the topological charge, which results in a mixing with the singlet quantum 
numbers as displayed in ( |164| ). 

10.5 Appendix E : Extended Bosonization. 

The use of the mean-field equation ( [I09| ) in Appendix A has allowed for a bosonization 
scheme that is trouble free. Indeed, if we were to carry a gaussian analysis around all the 
fields including the auxiliary field vr^ (and hence expand (|109|) , then instabilities show up 
along the scalar (vr^) and pseudoscalar (vr^s) directions. This, however, can be easily fixed 
through a generalization of (|13D to include fluctuations around the instanton densities. 
We start with a modification in (ITSf) 



1 = Jv7r^VP^exp(TTf Jdz P^{z)(n^{z) -^^^^^e^{z)^^ (165) 



that is inserted in the partition function (|^), where the term in the exponent 

n*/2 n* 
clearly couples scalar (pseudoscalar) glueballs to the quarks with strength gg {gps) (see 



m) below). 

Following the steps of appendix C, the mean field equation in tt^ reads 



~^ " Aimm{z) l-n*Ti±{z)/{AimN±{z)) ^ ' 
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Along with the contribution from the measure nin^, n_) (|^) the bosonized effective action 
reads 

5eff= - N,TTlogS-'[P^,a,x] - 2jdz—^TTjmP^{z) 
+ /..n,(.)log(iV,!n-^-^) 

+ Nf jdz [n+{z)+n-{z)^ +Sg (168) 
where Sq is the gluonic contribution (|126|) . The inverse quark propagator S^^[P^,a, x] 



in ([168D is 

S-i[P=t,a,x] = + S-1[P+,P- 



\ k"^ J n* 



X (^l + 7r, + «757rp,)(^l-^)yil4 (169) 

and its second term clearly exhibits glueballs coupling to quarks. The first term S^^ [P^, P^ 
in ( |169|) is given in (|17D . Following the saddle point approximation used in appendix C, 
we obtain the effective action as follows 

S,sKps,(J,x] = - N,TT\ogS~\P) +n*v(Nf - Tif^ ~ -) 
- N.Jdk 'K{%k)A+{k,mf,mg)7c3j{-k) 
+ N,Jdk7rfjik)A_ik,mf,mg)7rf,{-k) 

+ 5efr[x',xM (170) 

In the last two terms of (|170| ), we have lumped terms in o"^, arfg, x^ and X'^ps where tt^ 
(vTps) is decomposed in the scalar (pseudoscalar) unphysical basis (in flavour space) vrf^ 
(vr^l). Only diagonal terms vr^-^ (tt^/) mix with scalar (pseudoscalar) glueballs according 
to 



^eff[a^a7r,] = + j dk a{k) gli^k) a{-k) 

+ dka{k)J2(^-9s^^+ik,mf,mf))nfJ{-k) (171) 



f 
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where 

9:ik) = ^ + E (l - m„ m,)) - ^{g^ - 1)' (172) 

The point made at the start of this Appendix can now be appreciated. We see that if we 
set Qs = the resuh differs from what we had in Appendix C by the last term. The latter 
with its negative sign can cause an instability in the scalar glueball fluctuations. This is 
easily tamed by the use of the bosonization scheme discussed in Appendix C. 
For the pseudoscalar part, we have 

'S'e^r[x^X7^p.] = + J dk x{k) g^^{k) x{-k) 

+ fdkxik) J2(^-9ps'^A.{k,mf,mf)) t7rlf{-k) (173) 



where 

*w-2^:^l:^(i-||^-(*.'"/."/)) - w 

Again, if gps=0 the latter term may cause the fluctuations in the pseudoscalar glueball 
direction to be unstable. This is easily tamed by the bosonization scheme discussed in 
Appendix C. 

As a check, we clearly see that in the absence of fermion we recover the scalar (pseu- 
doscalar) gluonic contribution Sg for n± = (n + a ± x)/2. The expression (|171|) in terms 
of physical scalar meson field in Aq and Ag channels (flavour space) reads 

SesW^air,] = + Jdka{k)g^{k)a{-k) 

+ Jdk a{k) ^f{gl{k) n,o{-k) + gl{k) 7:,,{-k)) (175) 



and similarly (|173|) for the pseudoscalar sector reads 



'5e^r[x^xM = + dk x{k) gf,s{k) x{-k) 



+ jdk x{k) [glXk) Voi-k) + g^ik) r^s{-k)) (176) 

We have defined above for the scalar (pseudoscalar) part of the action 



9s,ps{k) = 1 - 9s,ps^^^-j^^ [ 2A±(fc, m, m) + A±(fc, m^, ) (177) 
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alpsi^^) = 9s,ps^j^'/^(^^±{k, rus, rus) - A±(/c, m, m)^ (178) 



The expression of and the mass gap equation of appendix C remain unchanged. 

• Scalar gluonic correlator Cff{x, y) 

The form of the scalar gluonic correlator is unchanged up to a constant term 

CM^ -y)^ fdk e^'^^-y^ f ^ + 4r E f iT-TT^ V -^)Y' (179) 

• Pseudoscalar gluonic correlator CpF^x^y) 
Similarly, for the pseudoscalar gluonic correlator 

1 1 / 1 2N.\\-^ 



• Pseudoscalar form factor 
For the pseudoscalar case, we have 

32n^ /I 1 ^ / 1 2A^,\ , , , 

^.(0) = TTiTT^r^TrTrTiT - + ^ E hr^TT-—- V - (181) 



2NcA_{0,m,m) \x* 2Nc ^ \A^{0,mf,mf) n 
as indicated in section 8. 

• Scalar form factor 

For the scalar case, we have 

^^Q, 327r^ / n 1 / 1 2NAy' 

' 2N^A+{0,m,m) \n^a^ 2iVc y V A+(0, m/, m/) J J ^ ^ 

as indicated in section 8. 

• Pseudoscalar coupling constant Qps 

In order to determine the strength gps{k — 0) of the pseudoscalar glueballs coupling 
constant to the quarks, we fit the experimental value of the mixing angle 9 — —20° when 
diagonalizing. The result gives an estimate for gpg {k—0 is understood) according to 

2 cot 26 {^9U^ - A.) + i) = y^iA - A.) + ^(^^ - i^) (183) 

We obtain gps = —7.025. 
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10.6 Appendix F : Unconnected Correlators 

We tabulate below the expressions for the unconnected correlator C}^{p) in the various 
channels 



7 






1 


2 ((1 - AO (1 - AO H - B,B,) 


A+{p) 


75 


-2 ((1 - AO (1- AO 11 + 5152) 


A-(P) 


7m 


2 (1 - AO 52 + ^ (1 - A2) 5i) 


A+{p) 


757/. 


2 (1 - AO 52 - ^ (1 - A2) 5i) 


A-(P) 




(p) = 






10.7 Appendix G : Expansion 

In this Appendix, we give the details leading to (jS^). Inserting the mass gap equation 
for small quark masses mi and m2 in the first term of A±(p, m, m) with the following 
approximations 

- mM M 

A ~ 5 ~ fl84) 

A;2 + M2-2mM' F + - 2mM ^ ' 

we obtain (denoting A±{p) = A±{p,m,m)) 

n , . f d^k (fciMa ± A;2M0 



n f 
A±{p) = ^ {miXi + 7712X2) + J 



{27:Y {kl + Ml) {kl + Ml) 



- (mi + 7712)/ 7 — -7—^ — T- — 185) 

At small momentum, and using (A;iM2 - A;2Mi)^ = p^Ml + {k.p)"^ (Mf - 2MkM'^/k), we 
obtain 
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2 f d'k {Ml - kM,M'j2 + Pil4V4) 
^^^^ " ^ A2vr)^ (A:^ + M|)^ 

(186) 

To illustrate the fact that the approximations used for A and B fail in the case of a strange 
quark mass, we show the plots of A (Fig. 23) and B (Fig. 24) for both the up (a) and 
strange quark (b). The solid line is the unexpanded result, while the dotted line is the 
expanded one. 



10.8 Appendix H. Outline of the Numerics 

In this appendix, we sketch how the numerical calculations were performed. First, we 
solve the integral equation ( pUD for A(m). With the constituent mass Mfc(m) fully known 
(p]), the propagator S{x,m) follows. The decomposition (^if) of S{x,m) lends itself to a 
straightforward numerical integration of Si{x, m) (^,^). The singular behavior at x = 
is contained in So{x,m). In p-space, each correlator is the sum of an connected part (|37|) 
and a unconnected part (|38|). To speed up the convergence of the numerical integration of 
(0), the free bubble diagram is removed by hand and later added. The evaluation of (|38|) 
is achieved in stages (numerator and denominator). Because we will later on numerically 
Fourier transform the p-space version of (0), great care is taken at low momenta since 
the reading of a meson mass is done at large distance. In x-space, the connected part of 
each correlator is directly evaluated from the x-space version of the propagator 5'(x,m) 
as 

C°{x) = -NcTr{S{x, m)-fS{-x, m)7) (187) 

For the unconnected correlator, we numerically Fourier transform its p-space version. In 
the pion case, we sum the two parts of the correlator and read the pion mass from the 
large distance behavior of x^^'^ times the correlator. The kaon unconnected correlator is 
exponentially damped by a factor of — with its pion analog and turns out to be 
of about the same order as its connected part. Therefore, we single out the unconnected 
part to read the kaon mass. 
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Figure captions 



Fig. 1. The constituent quark mass Mk{m) versus z = kp/2 {k being the momentum 
and p the average size of the pseudoparticle) for current masses m=0 (dashed curve), 
m=5 MeV (dotted curve) and m=10 MeV (solid curve), respectively. 

Figs. 2. The chirality flip (2a) and non-flip (2b) parts of the quark propagator 
(normalized to the free massless quark propagator) versus x (fm) for current masses of 
5 MeV (lower curve) and 10 MeV (upper curve), respectively. The squares are results 
of simulations carried out in for 128 instantons and 128 antiinstantons in a periodic 
box. The large distance behaviours are also shown for the chirality flip (2c) and chirality 
non-fliAp parts (2d). 

Figs. 3. The connected (3a) and unconnected (3b) parts of the correlator with 
arbitrary quantum numbers. The insertions correspond to the external background fleld 

as discussed in the text. The resulting connected (3c) and unconnected (3d) parts to 
leading order in l/N^ counting. 

Figs. 4. The connected part C^(p) of the correlator (normalized to the free and 
massless correlator). (4a) for the up (down) quark (m^ = nid = 10 MeV), and (4b) for 
the strange quark {ms=140 MeV) versus p (fm^-*^). The channels shown are scalar (S), 
pseudoscalar (P), vector (V), tensor (T) and axial- vector (A), respectively. 

Figs. 5. The same as Figs. 4 but for the unconnected part C}^{p) of the correlator. 
The tensor channel vanishes identically for the up and strange quarks as does the vector 
channel for the up quark. 

Fig. 6. The pion-correlator (normalized to the free and massless correlator) versus 
X (fm) for rriu = rrid = 5 MeV (upper curve) and rriu = rrid = 10 MeV (lower curve). 
The squares are the results obtained in using 128 instantons and 128 antiinstantons 
in a periodic box. The circles are the results obtained in 0] from cooled lattice gauge 
calculations. 

Fig. 7. The kaon-correlator (normalized to the free and massless correlator) versus 
X (fm) for rriu = = 5 MeV (upper curve) and for rriu = = 10 MeV (lower curve). 
The squares are the results obtained in using 128 instantons and 128 antiinstantons 
in a periodic box. 
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Fig. 8. In the pion channel, the connected C°(x, m„), unconnected C^{x,mu) and 
full correlators for 771^=10 MeV versus x (fm) are plotted in dashed, dotted and solid 
lines, respectively. 

Figs. 9. Large distance behaviour of the connected and normalized correlator 
C^{x,m). (9a) for up (down) and (9b) for strange quarks. 

Figs. 10. Large distance behaviour of the connected correlator times a;^. (10a) for 
up (down) quark, and (10b) for strange quark, respectively. 

Fig. 11. In the pion channel, the total correlator times x^^^ versus x (fm) is plotted 
for mu=5 MeV (upper curve) and mu=10 MeV (lower curve), respectively. 

Fig. 12. In the kaon channel, the connected correlator times x^^^ versus x (fm) is 
plotted for ms=UO MeV. 

Fig. 13. The connected correlator (normalized to the free and massless correlator) 
for the rj (00,88,08) versus x (fm). 

Figs. 14. The unconnected correlator times x^^^ for the (00) channel (14a) and the 
(88) channel (14b) versus x (fm) is plotted with (upper curve) and without (lower curve) 
a topological susceptibility. 

Fig. 15. The scalar connected and unconnected sigma meson-correlator (normalized 
to the free and massless correlator) versus x (fm) for rriu = fnd= 10 MeV are respectively 
plotted in short dashed and long dashed lines. The solid line represents their sum. The 
squares are the results obtained in using 128 instantons and 128 antiinstantons in 
a periodic box. The circles are the results obtained in 0] from cooled lattice gauge 
calculations. 

Fig. 16. The total correlator (normalized to the free and massless correlator) in the 
p meson channel versus x (fm), for = = 10 MeV. The squares are the results of 
36[| using 128 instantons and 128 antiinstantons in a periodic box. 

Fig. 17. The total correlator (normalized to the free and massless correlator) in the 
meson channel versus x (fm), for = m^; = 10 MeV. The squares are the results of 



36[| using 128 instantons and 128 antiinstantons in a periodic box. 

Fig. 18. The total correlator (normalized to the free and massless correlator) in the 
K* meson channel versus x (fm), for = = 5 MeV (upper curve) and = = 10 
MeV (lower curve). The squares are the results of [^] using 128 instantons and 128 
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antiinstantons in a periodic box. 

Fig. 19. The total correlator (normalized to the free and massless correlator) in the 
Ai meson channel versus x (fm), for m„ = rrid = 10 MeV. The squares are the results of 



3q| using 128 instantons and 128 antiinstantons in a periodic box. 

Fig. 20. The total correlator (normalized to the free and massless correlator) for the 
Ki meson channel versus x (fm) for rriu = 171^= 10 MeV. 

Fig. 21. The total correlator (normalized to the free and massless correlator) in the 
nucleon channel versus x (fm) is plotted for = = 5 MeV (upper curve) and for 
rriu = nid = 10 MeV (lower curve), respectively. The squares are the results obtained in 
using 128 instantons and 128 antiinstantons in a periodic box. The circles are the 



results obtained in |^ from cooled lattice gauge calculations. 

Fig. 22. The total correlator (normalized to the free and massless correlator) in 
the A channel versus x (fm) is plotted for m„ = m^^ = 5 MeV (upper curve) and for 
rriu = nid = 10 MeV (lower curve). The squares are the results obtained in using 128 
instantons and 128 antiinstantons in a periodic box. The circles are the results obtained 
in m from cooled lattice gauge calculations. 

Figs. 23. The coefficient A(k,m) versus z = kp/2 for up (23a) and for strange (23b) 
quarks. The solid line is the unexpanded result and the dotted line is the expanded result. 

Figs. 24. The coefficient B(k,m) versus z = kp/2 for up (24a) and for strange (24b) 
quarks. The solid line is the unexpanded result and the dotted line is the expanded result. 

Figs. 25. The scalar (pseudoscalar) gluonic correlator is plotted in Fig. 25a (b). The 
points are results from simulations in ([0). 

Figs. 26. In (a) and (b), we display the insertion mechanism involved in evaluating 
the gluonic form factor of a constituent quark. Figs (c) and (d) show the mixing that 
enters in the connected part. 

Figs. 27. The insertion mechanism for the fermionic form factor of a constituent 
quark. 

Figs. 28. The leading contributions to the fermionic form factor. 
Figs 29. The form factor in terms of loffe's current. 
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